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We search for strange pentaquark states that have been previously reported by other experiments 
- the 6(1540) + , S 5 (1860) , and ,~5(1860) — in 123 fb" 1 of data recorded with the BABAR 
detector at the PEP-II e + e~ storage ring. We find no evidence for these states and set 95% confidence 
level upper limits on the number of 0(1540) + and .55(1860) pentaquarks produced per e + e~ 
annihilation into qq and per T(4S) decay. For qq events these limits are about eight and four times 
lower, respectively, than the rates measured for ordinary baryons of similar mass. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



Ten experimental groups have recently reported nar- 
row enhancements near 1540 MeV/c 2 in the invariant 
mass spectra for nK+ or pK° s flElESSESSIlGil- 
The minimal quark content of a state that decays 
strongly to nK + is dduus; therefore, these mass peaks 
have been interpreted as a possible pentaquark state, 
called 6(1540)+. A single experiment (NA49) has re- 
ported a narrow resonance near 1862 MeV/c 2 in the in- 
variant mass spectra for £~ir~ and £~7r+ 11]. The 
minimal quark content of the £~7r~ final state is dssdu. 
Therefore, the latter two mass peaks have also been inter- 
preted as possible pentaquark states, named £5(1860) 
and £5(1860)°, with the latter being a mixture of ussuu 
and ussdd. On the other hand, a number of experiments 
that observe large samples of strange baryons with mass 
similar to that of the 9(1540)+ (e.g., 4(1520) -> pK~) 
see no evidence for the 6(1540)+ [l2j; a number of ex- 
periments that observe large samples of the nonexotic 
£~ baryon do not observe the £5(1860)" or £5(1860)° 
states |l2j . 

We report the results of inclusive searches for 6+ — > 
pK®, —>■ 3~tt^ , and £°^£~7r+ in e+e~ annihila- 
tion data, where we expect equal production of the charge 
conjugate states; their inclusion is implied throughout 
this Letter. The data were recorded with the BABAR de- 
tector ll| at the PEP-II asymmetric-energy e+e~ stor- 
age ring located at the Stanford Linear Accelerator Cen- 
ter. The data sample represents an integrated luminosity 
of 123 fb -1 collected at an e+e - center-of-mass (CM) en- 
ergy at or just below the mass of the T(4S) resonance. 

The BABAR detector is described in detail in Ref. 
We use charged tracks reconstructed in the five-layer sil- 
icon vertex tracker and the 40-layer drift chamber. The 
charged-particle momentum resolution is {<t(pt) Ipt) 2 — 
(0.0013pr) 2 + 0.0045 2 , where pt is the momentum trans- 
verse to the beam axis measured in GeV/c. Particles are 
identified as pions, kaons, or protons with a combination 
of the energy-loss measured in the two tracking detectors 
and the Cherenkov angles measured in the detector of in- 
ternally reflected Cherenkov radiation. We use all events 
accepted by our trigger, which is more than 99% efficient 
for both e+e~ — > qq and e+e~ — > T(4S) events. 

To evaluate the efficiency and mass resolution for re- 
constructing pentaquarks, we simulate pentaquark sig- 
nals with the JETSET 14] Monte Carlo generator by 



substituting a particle with the mass, width, and decay 
mode of a hypothetical pentaquark for an existing baryon 
already simulated by JETSET. We use large control 
samples of known particles identified in data to correct 
small inaccuracies in the performance predicted by the 
GEANT-based 0] detector simulation. The invariant- 
mass resolution for the decay modes studied in this anal- 
ysis ranges from less than 2 MeV/c 2 to approximately 
8 MeV/c 2 , depending on the final state and the momen- 
tum of the pentaquark candidate. 

We reconstruct 6+ candidates in the pK® decay mode, 
where if ° — > 7r+7r~ . A sample of K® candidates is ob- 
tained from all pairs of oppositely charged tracks we iden- 
tify loosely as pions (with more than 99% efficiency and 
[70-90]% rejection of K and p depending on momentum) 
that pass within 6 mm of each other. Each candidate 
is required to have a reconstructed trajectory passing 
within 6 mm of the interaction point (IP) in the plane 
transverse to the beam direction and within 32 mm of the 
IP along the beam direction, a positive flight distance, 
defined as the projection on its momentum direction of 
a vector from its point of closest approach to the beam 
axis to its decay point, and an invariant mass within 
10 MeV/c 2 of the nominal K® mass. We also require 
the helicity angle Oh of the candidate, defined as the 
angle between the 7r+ and the 7r+7r~ flight directions in 
the 7r+7r~ rest frame, to satisfy \cos0h\ < 0.8, reducing 
background from A decays and photon conversions. 

We combine these K® candidates with tracks we iden- 
tify as p or p (with [55-99]% efficiency and [95-99]% rejec- 
tion of 7T and K) that extrapolate within 15 mm (10 cm) 
of the IP in the plane transverse to (along) the beam di- 
rection. The invariant-mass distribution of pK® pairs in 
data is shown in Fig.^ No enhancement is seen near the 
mass of the reported 6(1540)+ (inset in Fig. 0. There 
is a clear peak containing 98,000 entries at 2285 MeV/c 2 
from A+ — > pK®, with a mass resolution below 6 MeV/ c 2 . 

We consider several additional criteria that might re- 
duce background to a pentaquark signal. Increasing the 
required flight distance of the candidates increases 
the A+ signal-to-background ratio, but does not reveal 
any additional structure. We also tried requiring at least 
one K~ and/or p candidate in the event. The A+ signal 
is still visible and there is no sign of a pentaquark peak. 

To enhance our sensitivity to any production mecha- 
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FIG. 1: Distribution of the pK® invariant mass for combina- 
tions satisfying the criteria described in the text. The same 
data are plotted for the full kinematically allowed pK® mass 
range and, in the inset, with statistical uncertainties and a 
suppressed zero on the vertical scale, for the mass range in 
which the B(1540) + has been reported. 



nism that gives a pK° s momentum spectrum in the CM 
frame (p*) different from that of the background, we split 
the data into ten subsamples according to the value of p* 
for the pKg candidate. The ten p* ranges are 500 McV/c 
wide and cover values from to 5 GeV/c, the kinematic 
limit for a particle of mass 1700 MeV/c 2 . The back- 
ground is lower at high p*, so we are more sensitive to 
mechanisms that produce harder spectra. There is no 
evidence of a pentaquark signal in any p* range. 

We quantify these null results for a + mass of 
1540 MeV/c 2 . We fit a signal-plus-background function 
to the pK® invariant-mass distribution for candidates in 
each p* range. We use a p-wave Breit-Wigner lineshape 
convolved with a resolution function derived from the A£ 
data and simulation. The latter is a sum of two Gaussian 
distributions with a common center and an overall root- 
mean-squared-deviation (RMS) ranging from 2.5 MeV/c 2 
at low p* to 1.8 MeV/c 2 at high p*\ this is narrower than 
the resolution due to the proximity of 1540 MeV/c 2 
to pK° s threshold. The best upper limit of 8 MeV/c 2 Q 
on the natural width T of the + is larger than our pK® 
mass resolution, and T could be very small. Therefore, 
we use T = 1 MeV/c 2 and T = 8 MeV/c 2 in the fit and 
quote results for each assumed width. We account for 
broad structures (known and unknown resonances, reflec- 
tions) in the pK® mass distribution by using a wide mass 
range, from threshold to 1800 MeV/c 2 , and a seventh- 
order polynomial times a threshold function for the back- 
ground shape; seventh is the lowest order giving an ac- 
ceptable x 2 - 

For the nominal selection criteria, we find that in each 
p* range the fit quality is good and the signal is consistent 
with zero. We consider systematic effects in the fitting 



procedure by varying the signal and background func- 
tions and fit range; changes in the signal yield are negli- 
gible compared with the statistical uncertainties. Vary- 
ing the mass assumed for the + has effects consistent 
with expected statistical variations. The other selection 
criteria give similar results. Since the nominal selection 
results in the smallest absolute uncertainties after effi- 
ciency corrections, we use it to set upper limits on the 
production cross section. 

We convert the signal yield in each range of p* into a 
cross section by dividing by the reconstruction and selec- 
tion efficiency, the K® — > 7r + 7r~ branching fraction, the 
integrated luminosity, and the p* range. If the + de- 
cays strongly, we expect only two possible decay modes, 
nK + and pK°, with very similar Q values, so we assume 
B(Q + — > pKf) = 1/4. The efficiency for the simulated 
pentaquark signal varies from 13% at low p* to 22% at 
high p* . The efficiency calculation is verified by mea- 
suring the differential cross section for A£ production in 
the combination of qq (q — d, u 1 s, c) and T(4S) events 
represented in our data. 

The resulting differential cross sections are shown for 
T = 1 MeV/c 2 and for T = 8 MeV/c 2 in Fig. The 
error bars include the relative systematic uncertainties 
on the luminosity (1%) and efficiency (4.9% dominated 
by the uncertainties on track and displaced- vertex recon- 
struction efficiencies). We derive an upper limit on the 
+ production cross section for each p* range under the 
assumption that it cannot be negative: a Gaussian func- 
tion centered at the measured value with RMS equal to 
the total uncertainty is integrated from zero to infinity, 
and the point at which the integral reaches 95% of this 
total is taken as the limit. These 95% confidence level 
(CL) upper limits are also shown in Fig. [2J 
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FIG. 2: The measured differential production cross sections 
(symbols) and corresponding 95% CL upper limits (lines) for 
Q + (top) and (bottom), assuming natural widths of F = 
1 MeV/ c 2 (solid) and at the current experimental upper limit 
(open/dashed), as functions of CM momentum. 



We derive model-independent upper limits on the to- 



tal number of pentaquarks produced per qq event and 
per T(AS) decay by summing the differential cross section 
over the kinematically allowed p* range for qq events (en- 
tire p* range) and for B meson decays (p* < 2.5 GeV/c), 
respectively, taking into account the correlation in the 
systematic uncertainty. The central value and the 95% 
CL upper limit on the total + (plus O ) production 
cross section for the p* range from to 5 GeV/c are 
shown in Table [J Dividing this limit and the corre- 
sponding limit for the p* range from to 2.5 GeV/c by 
the cross section for e + e~ — ► qq and for e + e~ — > T(4S), 
respectively, we calculate limits on the number of pen- 
taquarks per event, given in Table D For the maximum 
width (r = 8 MeV/c 2 ), we obtain a 95% CL upper limit 
roughly a factor of eight below the typical values mea- 
sured for ordinary octet and decuplet baryons of the same 
mass 16]. 

We search, as well, for the reported £5(1860)"" and 
£■5(1860)° states decaying into a E~ and a charged pion, 
where E~ — » yl°7r _ and A — > pn~ . We reconstruct 
A° — > p-K~ candidates from all pairs of charged tracks 
that satisfy loose proton and pion indentification require- 
ments and pass within 6 mm of each other. The yl can- 
didate must have a positive flight distance from the IP 
and an invariant mass within 10 MeV/c 2 of the nomi- 
nal A mass. These A° candidates are combined with 
an additional negatively charged track passing loose pion 
identification requirements to form S~ candidates, which 
are required to form a good vertex, to have a positive 
flight distance from the IP, and to have an invariant mass 
within 20 MeV/c 2 of the nominal E~ mass. The flight 
distance of the A° candidate from the A°ir~ vertex is 
required to be positive. This selection yields 290,000 
E~ candidates with a peak signal-to-background ratio 
of 23 in the A°ir~ mass distribution. Finally, we com- 
bine the S~ candidates with an additional charged track 
consistent with coming from the IP and passing loose 
pion identification requirements. The cosine of the angle 
between the reconstructed E~ trajectory, extrapolated 
back to the IP, and the additional track is required to be 
less than 0.998. This last requirement is especially im- 
portant, since the S~ is charged and has a long lifetime; 
if it has a long flight distance, it can produce a recon- 
structed track that, if combined with itself, forms a false 
peak in the invariant-mass distribution. The reconstruc- 
tion efficiency for the simulated pentaquark signal varies 
from 6.5% at low p* to 12% at high p* . 

The invariant-mass distributions for S~tt~ and for 
S~ir + combinations are shown in Fig. EJ In the i;"7r + 
mass spectrum, we see clear peaks for the S"(1530)° 
and S' c (2470) baryons, but no other structure is visi- 
ble. There are no visible narrow structures in the ,5*"7r" 
mass spectrum. 

As in the + search, we divide the S~tt~ candidates 
into ten subsamples according to the p* value of the can- 
didate. We find no sign of a pentaquark signal for any 
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FIG. 3: S~ir (black) and 5~w~ (gray) invariant-mass dis- 
tributions. The same data are plotted for the full kinemati- 
cally allowed S _ 7r mass range and, in the inset, with statis- 
tical uncertainties and a suppressed zero on the vertical scale, 
for the mass range in which the £5(1860)"" and S5(1860)° 
have been reported. 



range of p* . We fit a signal-plus-background function 
to the £~tt~ invariant mass distribution, for each p* 
range. Here no broad resonances or reflections are ev- 
ident, and we perform simpler fits over a E~Tt~ mass 
range from 1760 to 1960 MeV/c 2 using a linear back- 
ground function. The resolution function is derived from 
the £(1530)° and S" C (2470) Q signals in data and simu- 
lation, and is described by a Gaussian function with an 
RMS of 8 MeV/c 2 . For the Breit-Wigner width we con- 
sider two possibilities, 1 and 18 MeV/c 2 , corresponding to 
a very narrow state and the experimental upper limit on 
the width [Tlj[. respectively. We fix the .5*5"" mass 

to 1862 MeV/c 2 . In all ranges of p*, the signal is con- 
sistent with zero. Systematic uncertainties on the fitting 
procedure are again found to be negligible compared with 
the statistical uncertainties, and variations of the 
mass and selection criteria give consistent results. 

We convert the measured yields for the — > E~ir~ 
decays into cross sections as described above for the 9 + . 
The efficiency determined from simulation is verified by 
measuring the differential cross section for the observed 
i; (1530)° signal. The average relative systematic uncer- 
tainty on the efficiency is 6.2% with a slight p* depen- 
dence, and is larger than that for the pR® mode because 
there are two displaced vertices and more particles in the 
final state. We have used a E~ir~ branching fraction 
of one-half for purposes of calculating cross sections and 
limits, under the assumption that the two-body modes 
E~tt~ and £~K~ dominate and have similar branching 
fractions. 

The measured cross section and 95% CL upper limits 
for Sj (plus production are shown in Fig. and 

Table D For T = 18 MeV/c 2 , the limit on the total 
production rate per qq event is roughly a factor of four 
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TABLE I: The measured total production cross section and 95% CL upper limits (UL) on the cross section and yield per event 
for O(1540) + and £5(1860) pentaquark candidates. The natural widths V = 8(18) MeV/c 2 refer to the upper limits on the 
widths of the 9(1540)+ (H 5 (1860)""), used in the fits. 





total production 


UL on total 


UL on yield per 


UL on yield per 


particle 


cross section (fb) 


cross section (fb) 


qq event (10 _5 /event) 


T(4S) decay (10~ 5 / ever rt) 




r = 1 T = 8(18) 


T = 1 T = 8(18) 


r = i r = 8(i8) 


r = 1 r = 8(18) MeV/c 2 


e+ + e~ 


-19±93 7±183 


171 363 


5.0 11 


18 37 


s 8 -- + st + 


-53±25 -93± 38 


25 36 


0.74 1.1 


2.4 3.4 



below the typical values measured for ordinary octet and 
decuplet baryons of the same mass 0] . 

We perform a similar search for £5 — ►£' _ 7r + ) finding 
no signal in any p* bin. Since many decay modes are 
kinematically accessible to such a state with a mass of 
~1862 MeV/c 2 and the branching fraction is unknown 
a priori, we omit this state from Tableland express our 
upper limit on the total production of and S 5 per qq 
event as 0.8 x lQ- 5 /B(~l -> S~ir + ), at the 95% CL. 

In summary, we have performed a search for the re- 
ported pentaquark states 9(1540) + , ^(lSGO) -- , and 
£5 (I860) in e + e~ annihilations. We observe large sig- 
nals for known baryon states but no excess at the mea- 
sured mass values for the pentaquark states. 
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